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Abstract 

We re-examine the matter radii of diffuse halo nuclei, as deduced from reac- 
tion cross section measurements at high energy. Careful consideration is given 
to the intrinsic few-body structure of these projectiles and the adiabatic na- 
ture of the projectile-target interaction. Using n Li, n Be and 8 B as examples 
we show that data require significantly larger matter radii than previously 
reported. The revised value for n Li of 3.55 fm is consistent with three-body 
models with significant ls-intruder state components, which reproduce exper- 
imental 9 Li momentum distributions following n Li breakup, but were hitherto 
thought to be at variance with cross section data. 

PACS numbers: 24.10.Ht, 24.50.+g, 25.10,+s, 21.10.Gv, ll.80.Fv 
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Reaction cross section measurements at energies of several hundred MeV/nucleon have 
been used to study the radial extent of matter densities of short lived exotic nuclei produced 
by fragmentation Extensive tables of deduced radii are now available in the literature, 
e.g. ||. Glauber theoretical methods [|J[5| have been the basis for these assignments, and in 
particular the approximation HQ in which it is assumed that the projectile and target nuclei 
present static density distributions |J whose geometric overlap determines the reaction cross 
section. To high accuracy, the deduced rms radii are found to be essentially independent of 
the details of the projectile density distributions assumed, e.g. |||7]]. The accuracy of such 
deduced root mean square (rms) radii is of considerable importance since they are routinely 
used as empirical measures in constructing, constraining and assessing theoretical models of 
halo structures for use in the interpretation of data. 

At the heart of the static density model is the neglect of correlations between the pro- 
jectile (and target) constituents, each projectile nucleon being assumed to carry the same 
single particle density This assumption would appear to work well for spatially local- 
ized nuclei such as 12 C ||. For weakly bound systems such as halo nuclei, however, the 
intrinsic few-body nature or granularity of the projectiles imply strong spatial correlations 
between the valence nucleons and the more localized core. At incident energies of order 800 
MeV/nucleon one must also consider the relevant timescales for a significant motion of these 
valence particles inside the projectile and that for the passage of the same particle through 
the target interaction region. In breakup studies narrow momentum widths are associated 
with these valence particles which have characteristic kinetic energies of order 10-40 MeV 
within the projectile |J. For this reason reaction models p3| , |l"T| make an adiabatic approxi- 
mation, freezing the position coordinates of the few-body projectile constituents during the 
interaction. Physical observables are then obtained by suitably averaging the resulting po- 
sition dependent reaction amplitudes over the relevant position probability distributions of 
these constituents. 

This few-body picture suggests a quite different description of the projectile-target inter- 
action and formulation of the reaction cross section. Consider for example 11 Li as a pair of 
neutrons bound to a 9 Li core. For an impact parameter b of the n Li center of mass, Figure 
1, such that the projectile static density (shaded circle) overlaps the target, many spatial 
configurations of the constituent bodies will not overlap the target. The expectation is that 
the valence nucleon (large b) contribution to the reaction cross section will be reduced or, 
alternatively, that the collision will appear more transparent than otherwise expected. Nish- 
ioka and Johnson (TT2] investigated related adiabatic effects on light-ion composite projectile 
(d, t, 3 He and a) cross sections in the energy range 100 < E < 350 MeV/nucleon. The ef- 
fects were very significant for the extended deuteron but small for the a particle. Estimates 
of the accuracy of the static density model for n Li were considered earlier by Takigawa et 
al. for a simplified two-body (di-neutron) halo density and at lower energies, where the 
adiabatic limit is expected to be less reliable. They demonstrated clearly the convergence of 
the two-body and static density models in the limit of tight valence nucleon binding. They 
concluded that static density calculations would indeed overestimate reaction cross sections, 
however the model used was too crude to allow a quantitative discussion. The overestima- 
tion of cross sections in the static density model was also recognised previously by Chulkov 
et al. 0. 

In this Letter we examine the quantitative implications of the few-body adiabatic de- 
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scription for deduced matter radii of halo nuclei. We take as examples two-body, 8 B and 
n Be, and three-body, n Li, systems for which cross section data are available for each com- 
posite and core ( 7 Be, 10 Be, and 9 Li) nucleus on a 12 C target at 800 MeV/nucleon. Data are 
also available for the nucleon- 12 C system, so that all projectile constituent-target subsys- 
tems can be interrogated and compared with experiment. In common with the analysis of 
Takigawa et al. Jl]| we will apply the static density approximation to the spatially local- 
ized core-target and valence nucleon-target subsystems. Additionally, the adiabatic (frozen 
coordinate) treatment of these constituents allows us to study carefully the implications of 
a realistic treatment of the two- and three-body nature of the projectile wavefunctions on 
calculated cross sections. 

In Glauber theory 0] the reaction cross section for projectile P is 



POO 

a R (P)=2n dbb[l-T P (b)} 
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where Tp(b), the squared modulus of the Glauber S'-matrix, is the transparency of the 
collision at impact parameter b of the projectile center of mass (cm). In the static density 
limit 



T| D (6) = exp 
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where ct^n is the free nucleon-nucleon (NN) cross section, at the relevant energy, appropriate 



for the projectile and target Jl4| with densities pp and pr, and the 

p f\b) = r dz Pi {vwr^) 



(3) 



are the z-integrated densities or thickness functions. Here only the projectile ground state 
density enters the calculation and few-body correlations, the granular nature of the projec- 
tile, does not enter explicitly. 

In the few-body adiabatic limit, the transparency function is |]15||, 



T^(b) = \(<S>%\Sc(bc)S v (b v )\$ 



n\ |2 



(4) 



where |$q) is the wavefunction for the relative motion of the n-constituent bodies in the 
projectile ground state, the bra-ket denoting integration over these internal coordinates. For 
a two-body (one valence nucleon+core) projectile the core-target and valence nucleon-target 
*S-matrices, in the static density limit, are 



Sc(h 



v) 



t; 



1/2 



S v (b v ) = S 1 (b 1 ) = T^(h) 



1/2 



(5) 



with be the impact parameter of the core and T^ D the analogue of Eq. @ for the nucleon. 
For a three-body (two valence nucleon+core) system, then of course 



S v {b v ) = S 1 (b 1 )S 2 (b2) , 



(6) 



where the co-ordinates, in the plane perpendicular to the beam direction, are shown in 
Figure 2. Eqs. (|l|) through @ are calculated exactly in the following for realistic two- and 
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three-body wavefunctions |$q). The explicit forms of the three-body wavefunction for n Li 
are given in [11 . 

We apply the formalism above to calculate reaction cross sections in the static density and 
adiabatic limits for the one- and two-neutron halo nuclei n Be and n Li, and the one-proton 
halo nucleus candidate 8 B, all on a 12 C target at 800 MeV/nucleon. The choice of energy 
and target was dictated by our wish to connect cross sections for all binary sub-systems with 
experiment. 

For all three incident nuclei, the static density calculations of the projectile-target (7p D ), 
core-target (Sc) and valence particle-target (S v ) sub-systems use the prescription for 
(i = P, C, N) of Charagi and Gupta ]TJ]. A Gaussian matter distribution is assumed for 12 C 
in all cases with rms matter radius (r 2 )^ 2 = 2.32 fm ||. With these inputs, and assuming 
Gaussian matter distributions for the core nuclei with radii (r 2 )g /2 = 2.30 fm, (r 2 )^ 2 = 2.28 
fm and (r 2 )^ 2 = 2.31 fm, we calculate reaction cross sections for the core-target subsystems 
a R ( 9 U) = 796 (796 ± 6) mb, a R ( 10 Be) = 813 (813 ± 10) mb and a R ( 7 Be) = 738 (738 ± 9) 
mb. The empirical values, in parentheses, are taken from 0. The deduced core radii agree 
with those of |§ within error bars. The calculated nucleon- 12 C cross section at 800 MeV is 
a R (N) = 231 mb which also agrees with experiment [16] within quoted errors. Thus each 
projectile constituent-target input to the few-body calculations, the Sc and S v , is consistent 
with independent empirical data for that binary system. 

Part (a) of Figure 3 shows the results of static density and adiabatic calculations for the 
11 Li+ 12 C system for a number of theoretical three-body wavefunctions of n Li. We show 
the calculated cross sections versus the matter rms radius calculated from the wavefunc- 
tion models. The horizontal band shows the experimental interaction cross section datum 
a( 11 Li) = 1060 ±10 mb [17 and the vertical dashed line the previously quoted matter radius 
(r 2 ){{ 2 = 3.10 ± 0.17 fm §. 

The (upper) open symbols are the results of the static density model and the (lower) 
full symbols those of the adiabatic calculations for each wavefunction model. The reduction 
in the calculated cross sections, or increased transparency of the projectile in the latter 
case, is immediately evident. From left to right the diamond symbols correspond to the 
P0 through P4 intruder s-wave (Faddeev) model wavefunctions of Thompson and Zhukov 
||18|| , with increasing rms radius. The extreme right hand point is a continuation of these 
model wavefunctions (P5) with a ls-state scattering length of —44 fm and 80% (lsi/2) 2 
probability. The upright and inverted triangles are calculations using the L6A pairing model 
wavefunction |Tj|, which in the static density picture fits the published radial value, and the 
weak binding potential Os-wave intruder wavefunction (Gl of [|TBJ) inspired by the work of 



Johannsen, Jensen and Hansen f20|| . The straight lines through these model points are to 
guide the eye. 

The results of these calculations are indeed dramatic. Whereas static density calculations 
suggest a matter rms radius of order 3.1 fm, as reported previously, a correct treatment of 
the n Li three-body character now suggests the halo is very much more extended and that 
(r 2 )\{ 2 = 3.55 ± 0.10 fm, firmly in the middle of the range of values generated by intruder 
state models which successfully reproduce empirical breakup momentum distributions [[Tj| . 

Part (b) of Figure 3 shows the results of similar calculations but for the one-neutron 
halo system n Be. Again the horizontal band shows the experimental cross section datum 
a( 11 Be) = 942 ± 8 mb and the vertical dashed line the previously reported rms matter 
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radius (r 2 }{{ 2 = 2.71 ± 0.05 fm ||. The results are qualitatively very similar to those of the 
three-body n Li case. The angled dashed line shows the static density calculations and the 
angled solid line and full symbols the adiabatic model results. In this case these lines connect 
a large number of calculations using simple two-body (lsi/2) cluster wavefunctions for n Be 
using binding potentials with a range of geometries and with depth adjusted to reproduce the 
single neutron separation energy 0.503 MeV. The solid symbols are the results of adiabatic 
calculations for n Be wavefunctions |2T| which include the effects of core ( 10 Be) deformation 



and excitation. The wavefunction with rms radius 2.92 fm, whose calculated cross section 
lies within experimental error bars, best reproduces the excited state spectrum of 11 Be. 
These wavefunctions generate cross sections which follow precisely the trend of the inert 
core calculations and suggest a revised matter rms radius of (r 2 ){{ 2 = 2.90 ± 0.05 fm. 

Finally, in part (c) of Figure 3 we consider the one proton-halo nucleus candidate 8 B. 
The previously reported value of (r 2 )^ 2 = 2.39 ± 0.04 fm |J was very close to that for 
7 Be, (r 2 ) l j 2 = 2.33 ± 0.02 fm || suggesting, in spite of the very small proton separation 
energy (0.137 MeV) that the last proton had rather limited extension. The experimental 
cross section for 8 B has recently been revised to cr( 8 B) = 798 ± 6 mb [^2| and is shown by 
the horizontal band on the figure. Using the static density model and a Gaussian density, 
in the manner of 0, we obtain a revised static density estimate of (r 2 ) 1 ^ 2 = 2.42 ± 0.03 
fm, shown by the vertical dashed line. The angled dashed and solid lines are the results 
of static density and adiabatic model calculations for a large number of two-body (0p 3 / 2 ) 
cluster wavefunctions for 8 B based on Woods-Saxon potential geometries. The diamonds 
use wavefunctions based on the often used cosh form cluster model interaction |2~3l and lie 



on the same lines. Although the differences between the model calculations are smaller than 
in the neutron halo cases, they remain very significant and suggest the rms radius of 8 B 
should be revised to (r 2 )^ 2 = 2.50 ± 0.04 fm, indicating quite significant extension of the 
last proton distribution beyond that of the core. 

In summary, we have reanalyzed experimental data of reaction cross sections for n Be, 
n Li and 8 B projectiles on a 12 C target at 800 MeV/nucleon using an adiabatic treatment 
of the internal coordinates of the two- and three-body projectiles. We verify that all binary 
channel inputs to the adiabatic model are consistent with the available experimental data for 
these independent systems. The granular structure of the projectiles implied by realistic few- 
body wavefunctions is shown to reduce considerably the calculated reaction cross sections 
and increase significantly the values of matter rms radii deduced from data when compared 
to static density estimates. 

We deduce matter rms radii for 11 Li, n Be and 8 B of 3.55 ± 0.10 fm, 2.90 ± 0.05 fm and 
2.50 ± 0.04 fm, respectively, representing increases of 14.5%, 7% and 4.6% over previously 
tabulated values. Our revised radius for n Li is now consistent with theoretical three-body 
models with a significant ls-wave intruder state component, which reproduce breakup mo- 
mentum distributions, but were hitherto thought to be at variance with cross section data. 
Our revised radius for n Be is also consistent with two-body models which include core 
excitation and reorientation effects. 

The increased transparency of the few-body structures presented here is quite general, 
has implications for the deduced radii of all such exotic systems, and suggests that a careful 
re-examination of all such data is necessary. The particular importance of these effects in 
extended three-body halo systems is exciting. In the case of 11 Li we show this to be of 
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importance in elucidating their structure and in bringing consistency between calculations 
and data for both breakup momentum distributions and reaction cross sections. 
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FIGURES 



FIG. 1. Schematic representation of the static density (shaded circle) and few-body adiabatic 
(frozen coordinate) treatments of the three-body projectile (P)-target (T) collision at impact pa- 
rameter b. In the spatial configuration drawn the few-body projectile does not overlap the target. 

FIG. 2. Definition of position coordinates, in the plane perpendicular to the beam direction, in 
the case of a three-body (two valence nucleon+core) projectile. 

FIG. 3. Calculated static density and few-body adiabatic reaction cross sections at 800 
MeV/nucleon incident energy as a function of projectile rms matter radius, for a 12 C target. 
Parts (a), (b) and (c) of the figure are for n Li, n Be and 8 B projectiles, respectively. Details are 
given in the text. 
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